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Abstract 

Discussed is mechanics of objects with internal degrees of freedom in 
generally non-Euclidean spaces. Geometric peculiarities of the model are 
investigated detailly. Discussed are also possible mechanical applications, 
e.g., in dynamics of structured continua, defect theory and in other fields of 
mechanics of deformable bodies. Elaborated is a new method of analysis 
based on non-holonomic frames. We compare our results and methods 
with those of other authors working in nonlinear dynamics (many of them 
refer to our papers [2D], [2T], [35], [5D]). Simple examples of completely 
integerable models are presented. 

Keywords: afRne invariance, affinely-rigid bodies, collective modes, internal 
degrees of freedom, nonlinear elasticity, Riemannian manifolds. 

1 Introduction 

In our earlier papers [JS] , [SU] (and reference therein) we discussed systems with 
collective and internal degrees of freedom ruled by the afRne and linear groups, 
first of all the metric-preserving groups, i.e., the isometry group and its homo- 
geneous part, i.e., the rotation group (orthogonal group). Roughly speaking, we 
were dealing there with rigid bodies, i.e., gyroscopes, and affinely-rigid bodies, 
i.e., homogeneously deformable gyroscopes, both in the flat Euclidean space. 
Such objects are interesting from the point of view of purely rational analytical 
mechanics in itself, and besides, they occur in various quite practical problems 
P] . [55] . [iP] . [12, [13], [IS], [SO], [il]. In continuum mechanics they occur as 
models of internal degrees of freedom, i.e., microstructure models in mechan- 
ics of complex bodies. Classical examples are micropolar continua of brothers 
Cosserat and micromorphic continua of Eringen [M] , [15] , [16] , [17] , [19] , [35] , 
|36| . Roughly speaking, they are continua of infinitesimal gyroscopes or homo- 
geneously deformable gyroscopes. They describe some granular media; another 
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application is the long-wave continuum limit of dynamics of molecular crystals. 
Such continua do not consist any longer of the Newton type material points; 
instead, their elementary constituents are material points with extra attached 
orthonormal or general bases (respectively the micropolar and micromorphic 
continua). Let us mention there is also plenty of other microstructure models 
where at material points some other geometric objects are attached like e.g. 
liquid crystals which are, roughly speaking, continua of infinitesimal rods (in 
continuous limit of course) [TU] , [H] , [IH] , [S2 ■ To understand properly the dy- 
namics of generalized continua, one must start from the well-defined dynamics 
of extended afhnely- rigid bodies or their systems in a flat Euclidean space. 
Below we consider something more peculiar, namely the dynamics of very small, 
essentially infinitesimal rigid or affinely-rigid bodies in a curved Riemannian 
manifold, or even, more generally, in a manifold endowed only with affine con- 
nection. One must answer seriously the question concerning motivation. Such 
models are mathematically very nice and interesting on the level of purely ra- 
tional analytical mechanics. But it is clear that we are living in physical space 
which in a very good approximation is Euclidean (in the scale of our every-day 
life generally-realistic effects are negligible). Of course, one can answer imme- 
diately that this may be an intermediary step towards constructing relativistic 
mechanics of continua [S], [55], [IS], [51]. But topic, although interesting, is per- 
haps esoteric from the point of view of every-day praxis. Let us notice, however 
that there exist some problems where the motion of small objects in curved 
spaces becomes a practically viable topic. First of all let us observe that in 
a deformed elastic medium we are given two spatial metric tensors: the usual 
metric tensor gij of the physical space, which is in a sense the manifestation of 
the vacuum state of the gravitational field, and the Cauchy deformation tensor 

Cy, m, ED, 

In the above formula 77 denotes the metric tensor of the material space (the 
reference metric tensor) , and are Lagrange (reference) coordinates expressed 
as functions of the Eulcr (spatial) coordinates x'. Perhaps more naturally looks 
the expression for the reciprocal contravariant tensor. 

The point is that for small, concentrated objects in the deformed continuum, 
like, e.g., defects (dislocations, disclinations, vacancies, interstitials [3], [M]), it 
may be not gij but just Cy that is felt as something like metric tensor. Similar 
situations have place in solid state physics, where due to certain collective phe- 
nomena, it is used not the usual, geometric metric tensor, but rather the tensor 
of effective mass which in certain situations even need not be positively definite. 
It is reasonable to expect something similar here. Of course, the space endowed 
with Cij as the metric tensor is still flat due to the non-separability (compati- 
bility) conditions [19], [54], because the curvature tensor built of C does vanish. 
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But it is well-known that the point defects feel the body manifold as curved 
Riemannian space. Similar concepts appear in theory of residual stresses. 
Second, let us put our attention on two-dimensional surface phenomena. The 
surfaces of real bodies are curved submanifolds on dimension two. One can 
expect that for many surface phenomena it is just the two-dimensional induced 
metric that is relevant for the dynamics. For example, the microstructure of 
the bulk of the body may generate some effective microstructure on the bound- 
ary, and no doubt that for the corresponding surface phenomena the induced 
metric may be more relevant that the true metric tensor gij . Even without any 
microstructure, the induced two-dimensional metric seems to be an important 
factor, e.g., for the surface waves. Certainly it is so in dynamics of the boundary 
membranes (films) of biological cells. 

Let as also stress some ecological problems, like the motion of pollution regions 
(spots) on the oceanic surface after damages of tankers. This is just the clas- 
sical example of the motion of (relatively) small two-dimensional drop over the 
curved spherical surface. One can also think about the motion of continental 
plates and similar geophysical problems [20], [15], [12] • 

Of course, it is convenient to begin any analysis from the motion of bodies in 
constant curvature space, like the spherical or hyperbolic (Lobachevski) space, 
because one can expect rigorous analytical solutions, just in virtue of the exis- 
tence of high symmetries. There are some advantages of beginning the analysis 
from the academic case of n-dimensional spaces restricting only then the con- 
siderations to physical values of n. 

Within generally relativistic context the motion of extended structured bodies 
was studied by Mathisson, Weyssenhoff, Tulczyjew and Papapetru [32], [57] . 
[37] : let us stress the contribution of the first three Polish physicists. On the 
non-relativistic level the problems of rigid motion in curved spaces were dis- 
cussed by Yugoslavian mechanician Stojanowic [56]. We follow general ideals of 
Kiinzle [26], [45] for whom such an object was described by linear frame attached 
to the material point. This is a mathematical idealization of very small objects 
when the linear size of the body is small in comparison with the size at which the 
spatial curvature changes remarkably. What we do is somehow related to the 
pole-dipole approximation due to Mathisson, Weyssenhoff, Tulczyjew and Papa- 
petrou. It is reasonable to expect that performing some more detailled analysis 
one can obtain more adequate higher multipole description where equations 
of motion contain not only Riemann tensor but also its covariant derivatives. 
Therefore, the situation will be particularly simple in symmetric spaces (where 
the covariant derivative of the curvature tensor vanishes) [25j. The analysis 
probably would have to be lead with the help of normal coordinates and make 
use of certain asymptotic power series expansion. 

Before going any further, let us begin with explaining why we use here the ap- 
proach based on the concept of infinitesimal rigid or affinely rigid bodies, at 
least as a first approximation. 

For extended systems of material points in a flat space, any two admissible 
configurations of the usual rigid body (gyroscope) are related to each other by 
some isometry transformation. The isometry group of the flat Euclidean space 
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of dimension n is a ^n{n + 1)- dimensional Lie group. However, in a generic 
Riemann space the typical situation is that the isometry group has dimension 
zero and consists of merely identity transformations. In any case this dimension 
k always satisfies the inequality 

k<^n{n + l) (3) 

and the maximal possible value is attained in constant-curvature-spaces, includ- 
ing of course the flat spaces as a special example [25] . But even in a non-flat 
constant curvature space there are some problems, namely the centre of mass 
concept is not correctly defined [IS] , [5S] , and because of this there is no splitting 
into translational and internal motion. It is so when dealing with extended bod- 
ies. So there is only one way to escape the problem: to concentrate the attention 
on sufficiently small bodies and go to the limit with their size. Then we obtain 
just the non-Newtonian material point with something attached to it. In our 
models this "something" is an orthonormal basis (infinitesimal gyroscope) or 
just a general basis (infinitesimal affinely-rigid body, i.e., an infinitesimal homo- 
geneously deformable body). Then mathematically everything becomes correct, 
and physically one obtains a good approximate description of small metrically- 
rigid or affinely-rigid bodies. 

To be complete with all these arguments, let us remind the concepts of isome- 
trics in Riemann spaces and affine transformations in spaces endowed with afSne 
connection. Let (M, g) be an n-dimensional Riemannian space; M is a differen- 
tial manifold and g the metric tensor on M. We say that the mapping ip of M 
onto M is an isometry if 

9 = ^*9- (4) 

Analytically, when ip is given by the dependence of the new point coordinates 
on the old ones a;*, this is given by 



Let now F be an affine connection, e.g., the Levi-Civita one built of 5, but 
not necessarily, it may be a general connection, even completely non-related 
to g. Let V denote the corresponding covariant differentiation, and Vx the 
corresponding directional covariant derivative along the vector field X. We say 
that the diffeomorphism ip of M onto M is an affine mapping if for any vector 
fields X,Y ox\ M the following holds: 

V(^.x) {V*Y) = (VxF) , (6) 

where 1^9* is an abbreviation for the action of p) on vector fields (thus, the 
covariant differentiation is transparent with respect to the action of p). Just 
like isometrics, in a general manifold with connection, the existence of affine 
diffeomorphisms is rather exceptional. 
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2 Degrees of freedom, kinematics, phase space, 
symmetries 



From now on the physical space M is a differential manifold endowed with a 
Riemannian structure, i.e., with some positively definite metric tensor field g, 
or, in certain situations, even with the weaker structure given by some affine 
connection T. One can also consider the double structure when geometry of M 
is given by some pair [g, T) where some relationship between g and T may exist 
or not. As stated in Introduction, we replace extended bodies by structured 
material points with attached linear bases. When those bases are by definition 
g- orthonormal, wc deal with the infinitesimal gyroscope; if they are general, 
our system is an infinitesimal afhnely rigid body (homogeneous deformable gy- 
roscope). One can also consider intermediate situations when some weaker 
constrains are imposed onto infinitesimal affine motion (incompressibility etc.). 
The attached bases describe internal degrees of freedom and give the symbolic 
description of relative degrees of freedom after the limit transition, when the 
diameter of the body tends to zero. Let our object be instantaneously placed at 
the position x £ M and some liiiea^r frame 6 — (^i? • • ■ •, ■ ■ ■ ; 

) be attached 

at that point, thus 

CAeT^M, A=I,...,n. (7) 

For infinitesimal affinely rigid bodies those are general linear frames, for in- 
finitesimal gyroscopes they are orthonormal with respect to the metric g, 

(cAles) = 5 i^A, es) = gi^e"^ ag^ b = rjAB, (8) 

where 77 is a fixed reference metric in the micromaterial space N . Usually, but 
not always and not necessarily we put = i?" and take to be the Kronecker 
delta, 

Vab = 5ab- (9) 

The small latin indices refer to the physical space as seen, whereas the capital 
ones refer to the micromaterial space N . Injections of the body into M, or rather 
ino tangent spaces T^M are given by: 

y'(i,a)=x^(t) + ^VWa'', (10) 

where is the position of the material point, the quantities ip^ k describe the 
relative/internal motion, and are Lagrangian coordinates in the microma- 
terial space A^. More precisely, the formula (fTO|) is an approximation which is 
valid in virtue of the infinitesimal character of the body; the better valid, the 
smaller is the body. One can imagine it as formulated in terms of some normal 
coordinates in an e-order size neighborhood of the point a; G M, when e tends 
to zero. 
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The configuration space of an infinitesimal afRne body is given by the principal 
fibre bundle manifold of linear frames FM , or to be more precise, by the con- 
nected component of FM, to exclude the nonphysical singular configurations, 
when the system e = (e, . . . , e^i, . . . , e) fails to be linearly independent. Accord- 
ing to the ideas of differential geometry, interpretation of FM as a principal 
fibre bundle means that the group GL (n, R) acts freely and transitively on this 
manifold [25] . 

e = (e, . . . , e^, . . . e) ^ eL = (. . . , cbL^a, ■■■) (H) 

(do not confuse this action with the left hand side action of GL (T^M) on T^M 
and the manifold of bases in T^M). More precisely, the subgroup of positive- 
determinant matrices GL"*" (n, R) C GL (n, M) acts freely and transitively on the 
connected component of FM . 

Equivalently, instead FM one can use the manifold of co-frames F*M, i.e., 
frames of covariant vectors (linear functions on T^M at all possible x € M). 
There is an obvious natural diffeomorphism between FM and F*M given by 
the duality mapping 

FM3e={...,eA,...)^e={...,e^,...), (12) 

where 

e^{eB) = {e^,eB)=e^,e'B^S^B. (13) 
Here the group GL {n, R) acts on the left according to the rule: 

e=(...,e^,...)^eL=(...,L-iV,...). (14) 

If dimM = n then dim_FAf = dimF*M = n{n + 1) ; this is the number of 
degrees of freedom. 

Any system of local coordinates on M gives rise to the obvious coordi- 
nates (a;*,e*^) on FM or {x^,e^i) on. F*M; here e*^ are components of 
with respect to a;*. For simplicity we do not distinguish graphically between 
and their pull-backs o tt, x^ o tt* respectively to FM and F*M. Here 
TT : FM M, TT* : F* M M are natural projections which assign to the 
frames G F^M C FM, € F*M C F*M their attachment points x & M 
If M is an affine (flat) space with the linear space of translations V, then ob- 
viously FM, F*M are canonically diffeomorphic respectively with M x F{V), 
M X F{V*); F{V) and F{V*) denote here the manifolds of frames in V and 
V* (more rigorously - their connected components). This is the byproduct of 
the fact that the tangent and cotangent bundles TM, T*M are isomorphic re- 
spectively with the Cartesian products M x V, M x V* . Then the degrees of 
freedom split naturally into translational (M) and internal {V) ones. If M is a 
curved manifold, this is no longer the case; there is no canonical identification 
with Cartesian products, and often there is no identification at all because of 
topological obstacles. But of course, translational motion in M is always well 
defined due to the projection tt : FM M. If the total motion in FM (always 
well defined) is given by a curve g : M ^ FM (M is the time axis), then the 
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translational motion is described by the projected curve tt o g : R ^ M . Trans- 
lational velocity is also well defined and given at the time instant i e M by the 
tangent vector 

v{t) = (tt o q)' (t) - (Ttt o (t) e T^ie{t))M. (15) 

Unlike this, neither the internal motion nor internal velocity are well defined in a 
bare structure- less manifold. If our internal degrees of freedom are constrained 
to gyroscopic ones, the configuration space becomes {FM, g), the manifold of g- 
orthonormal frames. Then the number of degrees of freedom becomes reduced 
to 

dim(FAf, g) ^ -n{n + 1), (16) 

i.e., n translational degrees of freedom, parametrized by spatial coordinates 
and 

in(n-l) (17) 

internal gyroscopic deegres of freedom, in a sense like in the fiat Euclidean space. 
In the latter case the configuration space would be identifiable with 

MxFiV,g), (18) 

or rather with its connected components; F(V, g) denotes here the submanifold 
of g- orthonormal linear frames in V . Some important problem appears here, 
which is still important when the full system of affine degrees of freedom is 
concerned. Namely, in virtue of constrains ([5]) the quantities a;*,e*^ fail to be 
functionally independent and they are not generalized coordinates any longer. In 
a fiat space one can help with the problem of defining independent coordinates 
in the following way: one fixes some particular orthonormal reference frame 
E in M and then "parametrizes" F{V,g) via pT|) . using orthogonal matrices 
L g SOiy^g). Those in turn are parametrized in any of well- known ways, 
e.g., using canonical coordinates of the first kind, Euler angles etc. How to 
follow this pattern in a curved manifold? The only natural way is to introduce 
some non-holonomic orthonormal reference reference frame E in M, i.e., some 
auxiliary field 

E^{...,Ea,...) (19) 

on M, where 

{Ea\Eb) = g{EA, Eb) = g^JE'AE^B = Vab (20) 

and usually we put jjab = Sab, obviously when g (as assumed) is positively 
definite. If {M,g) is curved (Riemannian), then E must be non-holonomic. 
Let us remind that in differential geometry the non-holonomy object of a field 
of frames E is defined by any of the following equivalent formulas [64]: 

[Ea, Eb] - n^AsEc, dE^ = \n^BcE^ A E^ , (21) 
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where, as usual are elements of the co- frame E dual to E, "d" denotes the 
exterior differentiation, and for any vector fields X, Y their Lie bracket 
is analytically expressed by 

this vector field is correctly defined, i.e., it does not depend on the choice of 
coordinates [25]. The object O vanishes if and only if E is holonomic, i.e., 
consist of vector fields tangent to some coordinate lines. The corresponding 
tensor field on M, 

S[E] = ]-n^BcEA ®E^ ® E^, (23) 



i.e., analytically 



is the torsion of the teleparallelism connection Ftoi [E] built of E, 

^'jk = 2 (rtel[£^]'j7c - ^tc\[E]\-i^ , (25) 

The connection T[E] is analytically given by 

d 

Ttci[E]'jk = E'a-^E'^j. 

Geometrically it is uniquely defined by the demand that the vector fields Ea 
are all parallel with respect to the corresponding afline connection, 

\7Ea = 0, A^l,...,n. (26) 

When E is non- holonomic, the torsion of r[£'] does not vanish, but the curvature 
tensor Tl[E] of any r[£'] is always vanishing. Incidentally, let us mention that 
the teleparallelism connection is directly related to the non-holonomic represen- 
tation of any possible afhne connection. Namely, if T is some afhne connection 
on M, then its non -holonomic representation with respect to E is given by 
coefficients T^bc such that 

VcEb = T^bcEa; (27) 

here V denotes the covariant differentiation in the F-sense. One can easily show 
that 

r^jfc — ^tci[EYjk = E^ A^^ BC E^ j E'" k- (28) 

In particular, 

Ttoi[E]^BC = 0; (29) 

obviously, we mean here non-holonomic coefficients of Ftoi [E] with respect to E 
itself. 
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Let us go back to the main problem of our non-holonomic description. Assume 
that at some time instant t S M the structured material point is placed at 
x{t) G M and its internal degrees of freedom have the attitude 

e{t)^{...,eA{t),...)eF^^t)M; eA{t) e T,(t)M. (30) 

The reference frame just passed then is given by 

Ex{t) ^ {■ ■ ■ , EAx(t)i ■ ■ ■) ^ Fx{t)M;; E Ax(t) & Tx(t)M . (31) 

Being attached at the same point x{t) e M the vectors eA{t) may be expanded 
with respect to the frame Ex(t)i 

eA{t) = EBx(t)L^ A{t). (32) 

In this way, if E is globally defined and kept fixed, one can interpret the curve 

R 9 t -> L(t) e GL(n,M) (33) 

as a correct description of the internal part of motion, and the configuration 
space FM (its connected component) may be identified with the Cartesian 
product 

MxGL(n,M) (M X GL+(n,]R)) . (34) 

This representation becomes particularly important when gyroscopic constraints 
are imposed. Then obviously the matrices L{t) become orthogonal (?7-orthogonal) 
and the configuration space is identified with 

M X SO{n, M) [M x SO{i^)) . (35) 

Such description is very important technically, because, just like in the flat 
space model, the group SO{n,R) (in particular 5*0(2, R), 5*0(3,18) may be 
parametrized in a variety of standard ways, and in this way some good, non- 
redundant, generalized coordinates may be constructed. 

But even much more: Even without gyroscopic constraints, generalized coordi- 
nates (a;% c^a) which are well-defined then, are non-useful in realistic problems. 
Unlike this, the representation (|34p is well-suited to various viable dynamical 
models. The idea is to use the polar decomposition of L or its two -polar decom- 
position (singular value decomposition) [JH] , (JS] , [SO] ■ In resulting coordinates 
some physically viable dynamical models are analytically treatable. Because of 
this the description based on the use of auxiliary non-holonomic frames E is 
both technically useful and geometrically interesting. It enables one to reduce 
many curved -space- formulaes to ones formally very similar to expressions ap- 
pearing in mechanics of extended bodies in flat spaces [2D] , [31] . As mentioned, 
in a bare structureless manifold M there are two kinds of well-defined velocities: 
generalized velocity of the total motion in FM (or {FM,g)), p'it) parametrized 

by ( . . . , . . . , ^e*A, ■ • ■ ) and translational motion velocity v{t) — [tto p)'{t), 
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i.e., P3|) parametrized by ^. . . , . . The quantities ^e*yi fail to be well- 
defined internal velocities;moreover, they are not components of vectors in M. 
If some reference frame E is fixed, the time derivatives of L-matrices with com- 
ponents 

V(rl)^s(t):-|i-^BW (36) 

are well-defined measures of internal velocity. Usually it is more convenient to 
use Lie-algebraic objects denoted by firi, r^ri and given by formulas: 

f^rl^B = (l^^c) i-'^B, (37) 
^rl'^B = L^^^c-j^L^ B = L^^^C^yI" dL^ B- (38) 

Obviously, they are non-holonomic velocities, i.e., they are not time derivatives 
of any generalized coordinates. The label (rl) means relative; they describe 
internal motion with respect to the fixed reference frame E. Because of this 
essential dependence on the calculation-helping auxiliary quantity they are per- 
haps not very convincing, both physically and geometrically. Indeed E does 
not describe any real geometry. And it would be rather difficult to construct 
convincing models of kinetic energy (Riemannian geometry of the configuration 
space) using just these quantities and basing on the analogy with motion of 
extended bodies in flat space. Much more natural measure of internal motion 
velocities is based on covariant derivatives, because the affine connection jk is 
a "real" geometry of M, not an auxiliary analytical tool. Usually, although not 
necessarily, V^jk is the Levi-Civita connection built of the metric tensor gij 
on M: ^ 

r*ife = {]k} = {9m],k + gmk,j - gjk,m) ■ (39) 

Let us remind that is uniquely defined by the condition that the metric 
tensor is parallel, i.e., that its covariant derivative with respect to {j^,} vanishes, 

Vfcffy = 0, (40) 
and, in addition, it is symmetric, i.e., torsion-free: 

rv = rfc,- (41) 

In certain problems it may be useful to admit a weaker relationship between 
affine connection and metric structures, e.g., Riemann-Cartan space, Weyl space, 
Riemann-Cartan-Weyl space, or no relationship at all. There are also interest- 
ing models when only affine connection but no metric structure is assumed in 
M. 

Affine connection enables one to define vectors of internal velocities V"^ a'- 

V'a - ^/a - j/a + r,k {x{t)) e'Ait)^. (42) 
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This is simply the covariant differentiation of the attached vectors Ca along the 
curve of translational motion in M; the latter is described by the time depen- 
dence of spatial coordinates a;* {t) . When (M, F) is a flat affine manifold (vanish- 
ing curvature and torsion tensors), then obviously, FM trivializes to M x F{V) 
and the quantities -^e."^ a reduce in affine coordinates to usual derivatives -^e^A- 
In a general manifold with non-flat affine connection, the above quantities a 
(|42|) are non-holonomic velocities, i.e., they are not time derivatives of gener- 
alized coordinates. And the total covariant velocity (...,!)%...;..., V^a, ■ ■ ■) 
with components 

(.-.,»-..-;.-.i"..-.) = (.-..^-..-.^...-) («) 

is non holonomic in this sense. Let us now express these velocities in terms of 
the description based on the use of auxiliary (usually g-orthonormal) reference 
field of frames E. 

For any virtual motion, not necessarily one satisfying any equations of motion 
we have in virtue of (l32l) 



Va = ^^CA = {§iExWB^ L'^Ait) + E.,f^^)B^^L^A{t). (44) 

Obviously, from the point view of geometry of M, A{t) are scalar quantities, 
thus their covariant derivatives reduce to the usual ones, 

^ (45) 

The projected curve K 9 i ^ x{t) describes translational motion in M (analyt- 
ically, x'^{t)), and the first factor in the first term of will be expressed by 
the F-covariant derivative of the field E: 

■J^Ex(^t)B = {^zEB)x(^t) = \^cEB)x(^t) E "(*)"^ 

In analogy to the flat-space theory of extended bodies we shall use affine veloci- 
ties, i.e., what Eringen used to call gyration ^17) in this theory of micromorphic 
continua. 



^'j - ( §-/Aj e^, h^B - e\^/B. (47) 

Roughly speaking, this is respectively spatial and co-moving representation of 
the same quantity, 

n'j ^e^n-^Be^j. (48) 

If motion is (/-rigid, i.e., the frames ca are constrained to be orthonormal, these 
tensors are skew-symmetric with respect to the corresponding spatial and ma- 
terial metrics, 

= -n/ = -gjkg^'n^ (49) 
h^B = -^B^ - -riBcV^^'n^D (50) 
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and become the usual angular velocities. In exceptional but physical dimension 
n = 3 they are identified (as all skew-symmetric tensors) with axial vectors 
referred to as angular velocity vectors. 

We shall use the following notation for objects appearing in formulas (|44j|48p : 



a-i-^is = L-^^cj^L^B. (51) 
f^d/s = i-i^Fr^Dci'^Bi^G^^, (52) 
where are co-moving components of translational velocity vector, 

.« = e<=.f (53) 

The meaning of labels "rl", "dr" is respectively "relative" and "drive" .firi de- 
scribes the affine velocity with respect to the just passed reference frame i?a;(f) G 
Fx{t)M . r^dr is the affine velocity (angular velocity and strain rate) of the frame 
E itself as seen by the moving observer. This quantity describes the contribu- 
tion of translational velocity to the total b- This method of non-holonomic 
frames was intensively used by M. Zorawski in his book and papers on defect 
theory; |64| and references therein. 

Canonical momenta conjugate to generalized coordinates x*, (f^ a will be de- 
noted by Pi, p'^i- In many problems it is convenient to use non holonomic 
momenta conjugate to non-holonomic velocities fi^j, il^^s, v^. By analogy to 
mechanics of extended bodies in flat spaces we shall use the terms affine spin 
(hypermomentum) respectively in laboratory and co-mowing representation, de- 
noted respectively by S'j, Tj^ b- As mentioned, the covariant velocities of the 
total motion are given by 

(...."•........VV...) = (...,f........;^eV...); (54, 

if r is non-Euclidean, they are non-holonomic, 

V^A = v'A+T\jke'Av\ (55) 
v^A = pA. (56) 

Their conjugate non-holonomic momenta {pi,P^i) satisfy: 

P^ - P^ - e^AP^kT^.^, P^^ = p^^, (57) 

and the spin quantities are given by 

E^, = c'aP^j, ^^b = P^^e'B = e^AS^se^j (58) 

Let us observe the characteristic duality: is a well-defined vector in M, 
whereas v'^a = jt^^A fail to be so. Unlike this P^i = p^i are components 
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of well defined co-vectors in M, whereas Pi themselves do not represent co- 
vector components in M. The all dualities mentioned here and underlying the 
definitions of pi, P^i, T,^b are meant in the sense: 

P^V' + P^^V'a = P^V' + P^^V'A = P^v' + ^'j^h = PaV^ + B^^ A (59) 

where 

Pa = Pre\. (60) 

After some calculations one can show that the usual Poisson brackets for the 
phase space coordinates have the form 

{x\ x^} = {e\, e-'s} = 0, {x\ a} = {P^., - 0, (61) 

{P\x^} = 0, {e*^,P^,} = 5',5^^, {x\Pj}^6',. (62) 

There is nothing surprising in those brackets; they look exactly as ones with Pi 
replaced by pf, (remember that P^i — p^i] cf. ([57)1 . However, (pT|) ([S^ is not 
a complete system and the Poisson brackets like {Pi, Pj}, {Pi, P^j}, {Pi,e^A} 
are different than their holonomic analogues, like, e.g., 

{k,P,} = 0, {p„p^^} = 0, {p„e^A} = 0, (63) 

etc. Indeed, using the expressions one can easily show that 

{PijPj} = ^^iR^kij, 

{P^.P^,} = -P^fer^,, (64) 

where B}kij is the curvature tensor in the convention 

P kij — r f^j i r }^i j ~\- r a^r r ki (^5) 

and comma denotes the partial derivative with respect to the indicated coor- 
dinates a:^ Obviously, in a flat affine space, ((64|) reduce to ((63)) . The Poisson 
brackets (IM|) mean geometrically that. Pi are Hamiltonian generators of the 
parallel transport of state variables. In a flat space they are just Hamiltonian 
generators of spatial translations. It is easy to find another Poisson brackets, 
also convenient in applications, and at the same time admitting a clear geometric 
interpretation. For example, 

{j:\,j:',}^s%j:\~d\^'',. (66) 

Here we easily discover just the well-known structure constants of the linear 
group. And no wonder, S'^j, which are referred to as affine spin (hypermo- 
mentum) components, are Hamiltonian generators of linear transformations of 
GL{TxM) in FT^M (in the manifold of linear frames, i.e., internal degrees of 
freedom [T], [2], [5], [5D], gS], gS]- Let us remind, if we have the field of non- 
singular mixed tensors in M , 

Mex^L.e GL(T,M) C L(T,M) ~ TI{T^M), (67) 
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it acts on our configuration space, more exactly on the manifold of internal 
degrees of freedom, according to the rule 

F,M 9 e =(..., BA, ...) ^ (••■, o eA, ... ), (68) 

i.e., analytically: 

(..., x^ , e^A, ...) ^ (..., , V,{x)e=A. ■■■)■ (69) 

Obviously, the functions L^j above are components of the field L with respect 
to local coordinates a;*, 

L = L'j{x)-^(g)dx^. (70) 

Obviously, this action is essentially local in the sense of a;-dependence of L'^j on 
x''. In general there is no coordinate system in which they would be constant. 
And if accidentally they happen to be so in some coordinates, this is an artefact; 
in other coordinates they will depend on x'^ (unless L^j has the form X5''j , where 
A is constant). 

The quantities of affine spin, i.e., S^h are Hamiltonian generators of (|68)) in 
action on internal degrees of freedom. So if = exp(a(a:)) a being a mixed 
tensor field on M then ([68|) . ([69|) become 

F^M 3 e = {. . . ,eA, ■ ■ ■ ) {■ ■ ■ ,exp{a{x))eA, ...), (71) 

i.e., analytically, 

(...,x^ ..; ...,e\, ...) ^ (...,x^ (exp(a(x)))%(a:)e^"A, ...). (72) 

This means that the phase-space functions F depending only on coordinates 
but not on their conjugate momenta, i.e., F{x, e) suffer the transformation rule 

F^{m)F, (73) 

where the differential operator £[a\ has the following form in action on such 
functions: 

£[a]F = {a'j(x)E^„F} = (74) 
Another system of Poisson brackets, 

{p,,sM = E',r^,-i]^r',, (75) 

show that Pi are Hamiltonian generators of parallel transports along coordinate 
lines in M, just like T,^j (as we have just seen) are basic Hamilton generators of 
the system of groups GIj(TxM). One can easily show that (and this information 
is already contained in the above ones) for any function F depending only on the 
configuration, i.e., on the FAf-variables (but not on Pi and P^i) the following 
holds: 

{j:'^,F}^^E^,F = -e\-^F, (76) 
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where 

w^ = eK^^,E^^ = eK^. (77) 

For ((76)) . ((77l) no direct analogue holds for functions depending on canonical 
momenta (excepting the flat-space motion). 

Let us write down the co-moving representation of the above brackets, i.e., 
Poisson brackets for quantities Pa-, S'^b: 

\Pa.Pb\ = ^"^lB^kab - IPkS^ab. (78) 

{e^s,Pc} =-i5B5^c, (79) 
{s^b,S^d}-<5^bS-4z5-<J-^dS^b, (80) 

{i]^-,E^b}-0, (81) 

where the curvature and torsion symbols with capital indices denote the co- 
moving components of the corresponding tensors, so they depend not only on 
the spatial coordinates but also on the along- fibre coordinates a'- 

R'-kab = e^.i^^afce-'ife^e^B, (82) 

S'^AB = e'^S'.^e' Ae'^B. (83) 

Let us observe both similarities and differences in the structure of the co-moving 
system brackets and that one based on ((64l [66l [75| . Namely ([78l [79l [80|) is very 
similar to the basic commutation relations for the affine group GAj (n,K). In- 
deed, on the right- hand sides of ((751 [751 [50)1 we "almost" recognize the structure 
constants of GA/(n,K) ~ GL(n,M) x M" (x denotes the semidirect product 

of GL(n,M) and M"). The proviso "almost" concerns the first subsystem, i.e., 
([78]). Its right hand sides vanish only if the connection F is flat, i.e., if we 
consider the usual afhne-space-problem. Surprisingly enough, on the right-hand 
side the torsion tensor appears, unlike in This fact is both mechanically 

and geometrically interesting in itself, however for us it is not particularly in- 
teresting, because, as a rule, we shall work in a Riemann manifold, where F 
is the usual Levi-Civita connection, with vanishing torsion. Nevertheless, the 
Riemann-Cartan space, where the torsion is admitted, is also of some interest, 
first of all when the linear defects (dislocations) are admitted [M]. But the 
Poisson brackets (|79p . ([50]) do not depend on geometry of M and correspond 
exactly to the structure constants of GAj (n, R). It is seen that Pa are Hamil- 
tonian generators of parallel transports along the A-th "ribs" of A (roughly 
speaking). And similarly YAb are Hamiltonian generators of (fTTl [T4)l . There 
are however some structural differences between the role of E'j and S'^b in this 
sense. Namely, one can meaningfully say that L in (|lip is constant, it is just 
the system of scalars in M. 
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Of course, one can consider "micromaterial transformations" where L is lo- 
cal like in gauge theories, i.e., depends on coordinates . But nevertheless 
there exists a finite - dimensional subgroup, isomorphic with GL(n, R) itself 
and given by constant matrix elements L^b- Unlike this, as mentioned, the 
constancy of L^j in ([M1[7D|) is ill-defined. For the general, x-dependent matrices 
[L^b] = [(exp^)'^i3] we have, in a fuh analogy to ((TTIIT^ 

F^M 3 e = (..., CA, ...) ^ (..., CBLixfA, ■■■), (84) 

...,e\{x), ...) ^ (..., e's(exp£)^^, ...) . (85) 

In analogy to (|73p . any function of generalized coordinates F{x,e) suffers the 
transformation rule 

F ^ {exp^[£]) (86) 
where again the differential operator ^[£] acts on F as follows: 

mF = {^^a(2:)S^s, F} - £^a{x){^^b,F}. (87) 

And one can easily show that for any function F depending only on coordinates 
a;*, e*A we have 

r)F 

{S%,F} = -i?%F = -e's^, (88) 
{Pa, F} = ~HaF, {P,,F} - H,F, (89) 

where [30] 

H.^e^HA^^-T^,.e^B^^. (90) 

It is seen again that Pi are Hamiltonian generators of parallel transports along 
the axis of i-th coordinates and Pa are generators of parallel transports along 
the A-th legs of linear frames e. 

Let us remind the geometric meaning of the above quantities. In geometrical 
formulation of the theory of afFine connection the vector fields l are so-called 
fundamental fields of connection acting vertically along fibres in FM. Similarly, 
Hl are basic horizontal vector fields. Their structure; let us repeat it [30] : 

^"--^'-J^' ''^-^'^{ir.-r^eU^) (91) 

implies that their dual system of covector fields, i.e., differential one-forms uj^ , 
9^ is given by: 

= e^, {de^L + Vjke^Ldx'') , 0^ = e^.dx'. (92) 

As usual, by duality we mean the Kronecker- delta structure of mutual contrac- 
tions of basic covariant covectors fields with basic vectors, so: 

{lu''l,E^b) - S^'bS^l, {uj'^l.Ha) = 0, (93) 
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(0^,i?^B)=O, (94) 

Let us observe the characteristic geometric duahty. In (|9ip l do not depend 
on the affine connection F. They have to do only with the action of GL(ri,R) 
along fibres (they manipulate with internal degrees of freedom). Unlike this, 
by the very definition of being horizontal, Hk depend explicitly on the affine 
connection. And quite dually [i^^l] defining horizontal fields as its kernel, is 
explicitly built of F. But 6^ are also F-independent; they have to do only 
with the bundle structure of internal degrees of freedom. Incidentally 9^ do 
not depend on the choice of coordinates, although apparently they are defined 
with the use of coordinates a;*, but they, as a matter of fact, do not depend on 
the choice of x*. They may be defined without any use of M-coordinates at 
all, but here there is no place for discussing such purely geometric facts |3Qj . 
Let us observe that the above objects (pij) . in spite of their seemingly 
abstract geometric nature are nicely interpretable in purely mechanical terms. 
Namely, let : R — > FM be a curve describing the total motion of our object 
and let Q{t) £ Tgi^^^FM be its generalized velocity at the same time instant 
t S K. One can easily show that evaluating the differential forms w^i, 9^ on 
i.e., performing the total contractions of covectors lo^ l, 9^ with the vector 
Q{t), one obtains just the object of affine and translational velocities like (|49|) . 
(|50|) . ((53|) . Indeed, let us remind that the components of (|43| represent vectors 
in FM but not in M (except the components of translational velocity). But 
evaluating cu^ l, 9^ , we obtain well-defined co- moving components of the affine 
and translational velocity. Indeed, 

(iv'^L^g) =h''L, (^?^,^?) (95) 

Similarly, let us introduce instead ([?T1 [5^ the corresponding expressions: 
F J - e kF Le i-e Agjjj, 



9' = eVfi*^ = dx\ 
This is also system of dual bases (fields of frames in FM), 

{9\F^)^0, {9\Hk)^5\ 



(96) 



(97) 



The difference in comparison with ((9T1 [92)) is that ([97)1 are assigned to the 
particular choice of coordinates in M and the corresponding induced coordinates 
in FM, whereas (|91] l92|) are "objective", i.e., coordinate- independent. But 
nevertheless, in analogy to ([95]). 



= {^9\e)=v^ = ^. (98) 
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Roughly speaking (this is a "joke"), one can write: 

' dt ' dt' 

The doubled g-skewsymmetric part of E'j and ?7-skewsynimetric part of T,^b 
are the usual, i.e., metrical, spin quantities. They are given, by analogy to 
angular velocities, (|49l [50]) by expressions: 

S'^j := - E/ = E', - g^'gjil^'k, (99) 

V^B E^B - = E-^B - T^^^r/sijE^c. (100) 

More precisely, S*' J is the usual spin and V^b is what for objects in a flat space 
was called by Dyson "vorticity". S^j and V^b are respectively Hamiltonian 
generators of left acting (spacial) and right-acting (micromaterial) rotations 
preserving respectively 5(2:) and 77. Obviously, we mean the action on internal 
degrees of freedom, and by rotations we meant that in (|68p is gx- orthogonal, 
Lx e 0{TxM,gx), and in mi [H m L is 77-orthogonal, L e 0(R",7y), i.e., 
analytically, 

gxab^x'^'iLj' j = g^ij, VabL^kL^ M — Vkm- (101) 

It is important to stress that unlike in ([58)) . for deformative motion, V^b are 
not co-moving components of 5''^: 

S',^e\V^Be''f, (102) 

obviously, the same negative statement is true for extended bodies in a flat afhne 
space [45]. [49]. 

Let us stress a very important point: In flat afhne spaces there exist concepts 
like the total and orbital affine momentum (hypermomentum) and total afhne 
momentum with respect to some fixed spatial point. The same is true for the 
usual (metrical) spin or angular momentum. If are coordinates of the radius- 
vector of the centre of mass position x G M with respect to the mentioned fixed 
origin G M, 

ox = x'e, (103) 

(ei,...,e„ are basic vectors in the translation space V of A/), then the total 
affine momentum with respect to is analytically given by 

rj^x'p, + ^j, (104) 

where obviously, the first term is the "orbital" affine momentum. Similarly, for 
the total angular momentum, i.e., the doubled skew-symmetric part of J*j, we 
have the splitting into the "orbital" angular momentum and spin S^j, i.e., 
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angular momentum of the body with respect to the instantaneous position of 
the centre of mass in M: 

y,^U,+S^j. (105) 
As usual, the doubled skew-symmetric parts are meant in the convention: 

Z'j = J' J - Jj' = J' J - 9'^9jiJ\, 

L'j = x'Pj - XjP' = x'Pj - g'^Qjix^Pk, (106) 
o j — Zj j — Zjj — j ~ y 9jl^ k- 

In general, only the total quantities are well-balanced, or in special symmetric 
cases - just conserved. Nothing like the splittings (|104[ llOSp does exist in gen- 
eral curved manifolds. Moreover, it is only internal quantities, i.e., ones related 
to internal degrees of freedom that is well-defined. Neither orbital nor total 
quantities do exist at all. The reason is that in a general curved manifold there 
is no well-defined concept of the radius-vector. 

Obviously, Green and Cauchy deformation tensor, G[e] e M"* (g) M"*, C[e] e 
Tx*M (g) Tx*M (where e e F^M) are given by the usual formulas, thus, analyt- 
ically: 

G[e]AB=9'^j{x)e\e'B, C[e],j ^ r^ABe"^ ,. (107) 

These are standard formulas well-known from the theory of extended afhne bod- 
ies in a flat space; just the restriction of the general definition to homogeneous 
deformations. However, in both cases there is some delicate point the over- 
looking of which may lead to serious mistakes. The trap is hidden in the too 
automatic use of Schouten tensor notation. Namely, in certain formulas one uses 
contravariant tensors obtained by the metrical shift of indices, and following the 
Schouten formalism one writes them simply as: 

G[e]^^ = V^''v''''G[e]cD, C[eY' = 9"'9''C[eU, (108) 

According to the usual convention, when some metrics rj, g are fixed, the kernel 
symbols of tensors with rj— and g— shifted indices are identical with the primary 
ones and it is only the level of the indices that encodes the kind of tensor one 
is dealing with. The same convention is used for the inverse metric tensors, so, 
to avoid the crowd of characters one used the symbols 

VAB. V^"", 5y, 9^' (109) 

respectively for the covariant metric tensors and their contravariant reciprocals. 
The point is however that in certain formulas one uses contravariant tensors 
reciprocal to those with coordinates G[e]^B, C[e]ij. And denoting their compo- 
nents by G[e]^^ , C[e]*^ following (|109p would be completely misleading because 
of the confusion with (jlOSp . Therefore, for such an object one has to introduce a 
new kernel symbol, e.g., G[e], C[e] where the following coordinate-independent 
conditions are satisfied: 

, ^ j40 - ^ ik 

G[e] G[e]cB=S^B, C[e] C[e]kj^S'j. (110) 
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Those conditions are independent of 77,3, and except some special values of e, 
the inequalities hold: 

, - - a 

G[e] ^G[e]^^, C[e] ^ C[ef . (Ill) 

This is one of exceptional situations when the crowd of characters is unavoidable. 
When postulating and discussing dynamical models one must use transformation 
properties of all the above- introduced quantities under mappings ([55]) , ([55]) and 
([M]). pT|). It is easy to see that G[e] transforms under PT|) as follows: 

G[eL{x)]AB=G[e]cDL^Aix)L''B{x), eeF^M. (112) 

If for any x € M, L{x) is an orthogonal matrix, L{x) E 0(IR.",r/) C GL(n, R) 
then C is invariant: 

C[eL{x)] = C[e], e e F^M, (113) 

however, there is no particular well-defined rule when L(x) is a general element 
of GL(n,R). 

Transformation properties under ([55]) . are dual to the above ones. So for 
any mixed tensor field T on M given locally by T^j{x) {T{x) £ GL{TxM)), we 
have 

C[r(x)e]y = C[e]fc;T(x)^r(a;)'j, e e f^.M. (114) 

If for any x G M, T{x) is an isometry in the g{x)-sense, T{x) G 0(Ta;M, (/(X)) 
then G[e] is invariant: 

G[T(x)e] = G[e], e e i^,M, (115) 

and there is no well-defined rule expressing G\T{x)e] through G[e] alone if T{x) is 
not an isometry. Those are exactly the familiar rules well-known from mechanics 
of extended affine bodies in a flat space. Similarly, deformation invariants for 
internal degrees of freedom are given by the classical formulas. So, for example, 
we introduce the mixed Green tensors, 

G[e]^B ■■= V^^G[e]cB (116) 
and one of possible choice of basic invariants is the following 

i^,[e] :=Tr(G[e]") , a = 1, (117) 

One can easily show that 

j^a[e] :=Tr(a[e]-'^) , a = l,...,n, (118) 

where 

C[eYr-=9'''C[e]kj- (119) 
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Due to the Cayley-Hamilton theorem, quantities [e] constructed according to 
the rule (|117p but with other values of a € Z may be expressed as functions of 
the above ones. 

In spite of the above complete analogy with mechanics of extended afhne bodies 
in flat space, it must be stressed, however that in transformation rules of kine- 
matical quantities some essential changes appear in comparison with the flat 
space theory. And this is again a kind of "trap". Indeed, careful calculations 
show that ([551 affect velocities as follows: 

'n^j = r,i7'„,r-i''", + v^'^ (Vfer%„) T-i", (120) 

where the "primed" symbols denote the T-transformed quantities and Vy is 
the directional covariant derivative along the translational velocity V . This rule 
becomes identical with that for flat space only if T is covariantly constant, i.e., 

VfcT%„ = 0. (121) 

The dual canonical momenta transform as follows: 

'S*, = T\.S^„^-l™,, (122) 

^ B — ^ B- 

It is seen that the rule for internal quantities (spatial and co-moving compo- 
nents of the affine spin) is identical with that for flat spaces. Unlike this, the 
translational covariant momentum suffers an additive correction linear in VT. 
The micromaterial local transformations (IMl) act as follows: 



'V = V\ 

'^'j = ff,+e^B(iVfci-'^c)e^,V^', (123) 

11 B = C" D^^ B + 1^ Cl^ B,kV , 

'S^s - L-i^cS%L^s, (124) 

^ J — ^ J ' 

where, obviously, comma denotes the usual partial differentiation of scalar func- 
tions. It is seen again that the afhne spin transforms in a "proper" way, i.e., just 
like in a flat space. Translational covariant momentum and velocity variables 
transform "properly" when L is constant, i.e., when we deal with transforma- 
tions PT|) . 



21 



3 Equations of motion 



The first step is to derive equations of motion for non-dissipative, Lagrangian- 
Hamilton systems. Later on one introduces the general models by admitting 
some auxiliary generalized forces of non-Hamiltonian nature, responsible for 
non-conservative phenomena. The simplest way is to use the basic Poisson 
brackets introduced above. Kinetic energy of extended affine bodies in a flat 
space [48], [49] suggests us, by the simple analogy, to use the following formula 

T = Ttr + Tint = '^9^JV'V^ + 1 g,^- ^'^1/^ b J^^ , (125) 

where, let us remind, the positive constant m > is the mass of the body (in- 
ertia of translational motion), and J^^ are components of constant, symmetric 
[J^^ — J^^) and positively definite micromaterial inertial tensor (inertia of the 
rotational and deformative motion). But when dealing with internal degrees of 
freedom J is a primary quantity, no longer the second-order (quadrupole) mo- 
mentum of the mass distribution. Although such an interpretation may be 
admissible from the point of view of the mentioned limit transition with the 
size of the body, there are situations when such a procedure is essentially inade- 
quate, e.g., when dealing with such objects like gas bubbles in fluids, etc.. They 
have some inertial properties as elementary observation shows, but it is hardly 
expected that the internal inertia in such " exotic" situations may be " derivable" 
from something like 

J^^ = j a^a^dnia), (126) 

where are Lagrange coordinates and the positive measure n describes the 
mass distribution. It may be interesting and instructive to rewrite (|125p in some 
modified forms, e.g., 

T = Ttr + Tint = '^GabV^v'' + ^GKLn''Afl''BJ'''' (127) 

where are spatial, thus variable, components of the inertial tensor, 

JM^' ^^'^VsJ^^. (128) 

After performing the Legendre transformations for Lagrangians of the form 
L = T — il{x^,e^A) (no generalized, i.e., velocity-dependent potentials like, e.g., 
magnetic ones) and expressing generalized velocities by canonical momenta, we 
obtain the following expression for the kinetic energy (|125l) : 

X = Ttr + Tint = ^a^'P^Pj + \JABP^^P''o9'\ (129) 
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where pi, P"^i are canonical momenta conjugate respectively to x', e^A and J 
is the covariant inverse of J, 

J^^'JCB - S^B- (130) 

Warning: J must be not confused with the "covariant 77-shift of J", i.e., 

Jab^VacVbdJ^'''- (131) 

More precisely, (|129p is related to (|125p by the following explicit expression for 
the Legendre transformation: 

dL dT dx^ 

Energy is given by the usual formula 

8T 8T 

E = i'j-+e^Aj^-L = T + U{x\eJB) (134) 
and after Legendre transformation it becomes Hamiltonian 

H ^1 + H{x\e^B)- (135) 
In analogy to (I127p the kinetic Hamiltonian ()129p may be written in the form 

1^%r + %nt = ^G^''paPB + ]:JaB^^C^''d&° (136) 

2m 2 

This expression is built of geometric quantities-Hamiltonian generators of im- 
portant transformations acting in the configuration and phase spaces. More 
precisely, it is a quadratic form of generators. Coefficients of this quadratic 
form depend on the configuration variables a;', c'a, i-e., on the position of the 
body in M and internal parameters. 

The simplest way to derive equations of motion is not to derive them from 
Lagrange equations. This would be the terrible, very strenuous work where 
mistakes are simply unavoidable and the structure of resulting equations is non- 
readable. Rather, at least for non-dissipative systems, the best procedure is to 
use the above-quoted basic Poisson brackets and to write equations of motion 
in the form 

^^{F,H}, (137) 



where H is the Hamiltonian, e.g., (|135i I136P and F runs over some complete 
systems of functionally independent functions on the phase space (some non- 
traditional in the sense of Darboux coordination of our phase space). Later 
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on, one may perform the inverse Legendre transformation and go back to the 
usual second-order equations of motion equivalent to Lagrange equations of 
the second kind. For certain reasons it is instructive to begin with a rather 
academic situation when there is no a priori assumed relationship between afRne 
connection T and metric tensor g in our space manifold M. Then as a matter of 
fact we are given two afhne connections: T^jk and Levi-Civita connection {*^} 
induced by g. Obviously, as usual, the difference of two affine connections 

ii\c r%, - {1} (138) 

is a tensor field (once contravariant, twice covariant) although r°-bc, {°6c} sep- 
arately have not this property. 

Performing the above Poisson-bracket-procedure (jl37p of deriving equations of 
motion and inverting Legendre transformation, one obtains the following sug- 
gestive system of second-order differential equations written in a balance-like 
form: 



Dv' 
m—— 
Dt 



m^'jkv'v^ -f S"„i?"„%w^' (139) 



Dt Dt + ' 



= -eV5^"%^^ J^^'^ + (140) 
Dt'^ Dt Dt 

where the meaning of symbols is as follows: 

1. the shift of indices (up-down) is meant in the sense of metric 



2. 



N^' = N\g^\ N\ = -e^A-^ = ^E\^ (141) 



(cf GZH), 



3. 



= g'^'Fu, Fk = -i/fcil = - ^ + r^fee^B^ (142) 



(cf (HOD). 

Let us mention that i/fcil are components of what is known in mechanics of 
principal bundles of frames as a covariant exterior differential Dil of the scalar 
function il. Note that if it is not projectable from FM to M (i.e., if it depends 
not only on but also on e*^), then Dii ^ c?il unless the affine connection 
F is flat. Let us notice also that in a flat afflne space (Euclidean space) our 
equations of motion (|139p . (|140p reduce to the well-known equations of motion 
of extended affinely rigid bodies [35] j [IHj- And also in curved spaces the main 
ideas of interpretation are similar to the mentioned ones. On the left-hand side 
of (|139[) we have covariant acceleration multiplied by mass. Indeed, — ^ 
is just the usual velocity of translational motion. On the right-hand side we 
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have the usual force acting on the material point. Similarly, in (|140p N'^^ 
is just what was called affine momentum, or hyperforce, influencing directly 
the motion of internal degrees of freedom. We have started from variational 
principle, therefore, these dynamical quantities F*, TV*-' have the very peculiar 
potential structure being build of derivatives of il. Incidentally, let us observe 
a very interesting relationship: 



However, once derived in this way, equations of motion (|139p . (|140|) may be 
easily generalized by admitting dissipative phenomena. Simply, we may allow 
the force F'' and hyperforce N'^^ to be completed by some completely non- 
potential terms depending on velocities w*, V^a or equivalently w% Q'^j Such 
additional terms may describe viscous friction, both in translational and internal 
motion. In this sense, without the particular model (|141II142| equations (|139[) . 
(|140p are quite general. 

Let us now concentrate on peculiarities of motion in curved spaces. In a sense 
the right-hand side of (|139p may be interpreted as th total force F/^j affecting 
translational motion and responsible for its covariant acceleration. But except 
the external force F*, this total force F^^^. contains certain additional terms. In 
an obvious way those additional terms represent the geometric force in the sense 
of coupling between geometry of M and the particle degrees of freedom, both 
internal and translational one. Similarly, in (|140|) besides the usual "external" 
term N''^ , on the right hand side there is an additional geometric term describing 
the coupling of internal degress of freedom with spatial structure. Notice that 
even if external translational force vanishes, = 0, translational motion is 
not geodetic. This is just the consequence of the mentioned coupling between 
geometry of M and mechanical degrees of freedom of the particle. If we insisted 
on approximate description of a small extended body, this phenomenon would 
have to do with geodetic deviation. Equations of motion (|139p , (|140p simplify 
in a remarkable way when we go back to the model where affine connection T 
is implied by metrical structure g, or at least when there exists some kind of 
compatibility between F and g. 

The most natural model, applicable both in defect theory [64] and in General 
Relativity is that of Riemann-Cartan space where metric g is assumed to the 
parallel under affine connection F, 



Then A'^bc becomes the so-called contorsion tensor K°'bc and it is easily shown 
in differential geometry that [M] 



Fk 




(143) 



Vfc5y = 0. 



(144) 



^ be — ^ bc-r Obc + ^cb 



a 



(145) 



where S is the torsion tensor of F 




(146) 
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All tensor indices are manipulated in the sense of g; in particular, contorsion is 
g-skew-symmetric in the first pair of indices, 

K\c = -Kb\ = -gug^'K^jc- (147) 

In this way F is controlled by two independent and a priori arbitrary quantities: 
g and S. Equations of motion (|1391 1140p simplify then to the following form: 

m— = T.\R\',v^ + 2mv''v'Sb,' + F\ (148) 

A^e^ bJ^"" ^ N^' ■ (149) 

In Riemann-Cartan spaces the curvature tensor is g-skew-symmetric in the first 
pair of indices, 

R^aij = —Ra'ij — —gacg^'^R'^dij', (150) 

Because of this (jl50|) y-skewsymmetry it is only the g-skew-symmetric part of 
E that really enters the formula (|148p . But this skew-symmetric part is just the 
half of internal angular momentum (spin) in the spatial representation, 

i ~ S/') = i (E^, - .g^'^g.^E™,.) = ^S^, (151) 



cf. (HH). Therefore (fT48)) becomes: 

= = 2mv-v'Sab' + ^S\R\',v^ + r (152) 
and the similar balance-like form of (|149|) reads 



JabE^E'-J +iV'^ (153) 
where the e-dependent inertial tensor Jab is given by 

Jab = JKLc'^ae'-b, (154) 

or, equivalently, in reciprocal contravariant terms 

r'' = eVe^ J^^. (155) 

Equations of motion (|152[ I153p are represented here as a system of balance laws 
for the kinetic linear momentum and affine spin p^, E*-', where according to 
Legendre transformation 

f = mv\ E*j = g^a^'^bJ'", E^^' = n^bJ'"- (156) 

Let us stress that unlike in the usual formula 

P^^g.,^J^^ (157) 
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the " coefficients" j''* in tfie non-holonomic representation (|156p depend on the 
internal configuration e. So, they are non-constant and are state-dependent 
even in the flat-space theory. Formally the internal part of dynamics (|149p . 
(|153p looks like in mechanics of extended afline bodies in Euclidean space. The 
difference is only that the usual time derivative ^ is replaced by the covariant 
one The main novelty, as mentioned above, appears in the dynamics of 
translational motion (|148p , (|152[) . Because now it is not only so that the " flat" 
^ is replaced by Namely, in addition to the "external" translational 

force two "geometric" forces appear on the right-hand side of (|148p . (|152p . 
One of them describes the direct coupling between linear momentum (thus also 
translational velocity) and torsion tensor. The other one describes the direct 
coupling between internal angular momentum (spin) and the curvature tensor 
R. Those geometric forces have the " magnetic-like" structure in the sense that 
they are 5-orthogonal to translational velocity, so they do not do any work. The 
duality: translations-torsion, rotations-curvature is well known in differential 
geometry [53], and on the mechanical level-in defect theory (dislocations and 
interstitials or voids) [3], [64]. This is a nice picture, however, let us observe 
that (|1481 1152p may be also written in shorter form, where the translationtorsion 
term apparently disappears, 

f ^-l^-^^^i^V.V +FN (158) 

where ^ denotes the covariant differentiation in the g-Levi-Civita sense. Ob- 
viously, if there is no torsion then ^ = • 

The above equations of motion are structurally similar to generally-relativistic 
equations for the pole-dipole particle. The mentioned generally relativistic 
models were investigated by Mathisson, Weyssenhoft, Papapetrou, Tulczyjew, 
Kiinzle and manv others [32] . [37] , [57] . The characteristic couplings: translations- 
torsion, rotations-curvature have a very deep geometric background and in this 
way, "geometric forces" on the right-hand side of (|1481 1152]) . although derived 
from some dynamical Lagrangian postulates, are based on some kind of a priori. 
We concentrated here on the Riemann spaces, when Vg — and torsion tensor 
does vanish. Riemann-Cartan space is in a sense more natural, although there 
appear some ambiguities in the definition of angular velocity and afflne velocity 
when torsion is admitted. We do not deal here with more general relationships 
between F and g, although some of them should be mentioned as interesting in 
some possible future investigations. Riemann-Cartan- Weyl space is defined as 
such one in which parallel transport preserves angels between vectors but not 
necessarily their lengths (roughly speaking parallel transport acts on vectors as 
a conformal transformation). Then we have 

Vfcffij = -Qkgi], (159) 

where 

Qfc = g^'^'Qra (160) 
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is traditionally referred to as the Weyl vector {Qk is then the Weyl co- vector). 
Then ^38]) has the form 

Si'jk = S',k + S.k" + 5^-, + {5',Qk + 6\Qj - gjkQ') , (161) 

where, as usual, tensor indices are shifted in the g-sense. If F is symmetric (no 
torsion) then obviously: 

ii'jk = S'jQk + 5\Qj - gjkQ' (162) 

and one is dealing with what is usually called the Weyl space. Riemann-Cartan- 
Weyl spaces and Weyl spaces are interesting in themselves and useful in defect 
dynamics [64j . However, the problem of how to define angular and affine ve- 
locity becomes here essential. In any case in this paper we concentrate on 
Riemann-, or sometimes- Riemann-Cartan spaces. In mechanics of extended 
affine bodies in fiat spaces one often discusses the problem of affine dynamics 
with additional constraints [IH] , [IS] , [SH] ■ The most important example is that 
of gyroscopic motion, when the body is metrically-rigid. By analogy we can, 
and in many problems just should, discuss the motion of infinitesimal gyroscope 
in a curved manifold. According to gyroscopic constraints the frame e should 
be (7-orthonormal during the motion; analytically: 

VABe^^{t)e^J{t) ^ g,, {x{t)) . (163) 

This means that affine velocity is permanently (;-skew-symmetric: 

= -n/ = -g,fer!^,9™, (164) 

just angular velocity in spatial representation. The corresponding equations of 
motion may be obtained from (|148II149]) . i.e., from (|152[ I153[) . even without 
the use of variational principle, just basing on the d'Alambert principle of ideal 
constraints. The power of generalized forces F% iV*-' on virtual generalized 
velocities w*, ri*j is given by 

^jfi^Fy + N'jn^,. (165) 

According to the d'Alambert principle, equations of gyroscopically constrained 
motion are obtained in the following way: 

1. on the right-hand side of unconstrained equations one introduces addi- 
tional reaction forces -Fr, Nfi which maintain constraints, 

2. one substitutes formally constraints equations, 

3. one assumes that Ffl,, iVfl do not do any work on virtual displacements 
compatibile with constraints: 

= F^y + NR\n^, = (166) 

for any v and any fl satisfying p64p . 
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The result is that does vanish and Nji is g-symmetric, i.e., its f;-skew- 
symmetric part does vanish: 

Fr' = 0, Nr', - Nr/ = Nr', - 5,fcAffl^„5"' = (167) 

or briefly: 

Fr' = 0, Nr'^ - Nr^' = 0. (168) 

This means that the effective reactions-free system of equations of motion consist 
of (|148p (equivalently (|1521 1158P , the skew-symmetric part of p49p (equivalently 
one of (|153p ) and constraints equations (|163p . The mentioned skew-symmetric 
part of internal equation 

(^e\^e'B - e^A^e^B^ J^^ = N^' - N^' (169) 
has the following balance form given by the skew-symmetric part of (I153P 

= N'^ - N^' (170) 

Dt ^ ' 

where, obviously, j are spin components (|99i I15ip . The skew-symmetric in- 
ternal hyperforce 

= N'^] _ ]Sf3i (171) 

is just the usual torque, i.e., moment of forces acting on the body. In n dimen- 
sions it is just a skewsymmetric tensor; the peculiarity of the physical dimension 
n—i is that it may be identified in a known way with the axial vector 91*. 
Just like in mechanics of extended affine bodies in flat spaces one can consider 
also other physically interesting constraints like incompressible motion, rota- 
tion -less motion etc. The corresponding equations of motion consist of (|152p 
i.e. (|158p and respectively the trace-less part of symmetric part etc. of (|149p 
i.e. (|153p . and, of course, equations of constraints themselves. 
The above balance form of equations of motion is very instructive and reveals 
geometric foundations of the model, flrst of all its symmetry properties and the 
corresponding conservation laws. Another problem is how to solve equations of 
motion and determine the phase portrait, at least qualitatively. And here the 
problem of using the non-holonomic reference frame E is crucial. There are two 
main reasons for that: 

1. As already mentioned, when gyroscopic or other constraints are imposed, 
the quantities (a;*, c^a), or equivalently [x^.e^i) are no longer independent 
generalized coordinates. Then the best way to introduce reasonable and 
computationally effective generalized coordinates is just to introduce an 
auxiliary field of (co)frames E, use representation ([5^ . confine matrices 
L to be elements of S'0(n,M) (gyroscopic model) or some other subgroup 
G C GL(n.M), and parametrize G with the use of some natural coordinates 
e.g., canonical coordinates of the first or second kind, Eulcr angels, some 
byproducts etc. 
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2. Even if no additional constraints are imposed on affine motion of internal 
degrees of freedom, it is rather a rule than exception that the natural 
coordinates (x*,e'yi) on FM are inconvenient and non-effective in study 
of realistic problems like elastic vibrations and their coupling with rigid 
rotations. For example, one deals often with isotropic problems when the 
potentials energy il is built of deformation invariants. The only reasonable 
procedure is then to use ([32)1 and express L in terms of left or right polar 
decomposition or two-polar decomposition (singular value decomposition). 
Incidentally, this enables one to use description as similar as possible to 
the one of extended bodies in flat space. 

For any L € GL(n, M) we have two version of the polar decomposition, some- 
times referred to as the left or right one: 

l = os^j:o, i: = oso-^, (172) 

where O G 0{n,]S.) is orthogonal and S ~ S'^ , S = E"^ are symmetric and 
positively- definite. Performing diagonalization of S with the help of some 
orthogonal matrix V G 0(n,M), 

S = VDV-\ (173) 

D being diagonal and positive, and denoting: 

U = OVeO{n,R) (174) 

we obtain the singular value decomposition (two-polar decomposition): 

L^UDV-\ U, VeO{72,R), £>- diagonal. (175) 

It is well known that the polar decomposition is unique, but the singular value 
decomposition suffers some kind of multivaluedness which is essentially harm- 
less if properly treated [iH] , 03] . 

In (|172p internal degrees of freedom are represented as consisting of two sub- 
systems: rigid body in n-dimensions and deformations; respectively ^n{n — 1) 
and ^n{n + 1) degrees of freedom. There are two possible representations of 
deformative modes, as seen in (|172|) . The symmetric-deformative objects have 
respectively to do with the Green and Cauchy deformation tensors. Namely, 
matrices of those tensors are given by 

G = L^L = S^, C = E-^ C-1=S2. (176) 

In the two-polar decomposition L consists of two fictitious rigid bodies repre- 
sented respectively by U and V; every with ^n{n — 1) degrees of freedom, and 
of n purely scalar deformations (stretching) which tell us only how the body 
is stretched, but without any information about orientations of this stretching 
both in physical space and in material of the body. The quantities U, V describe 
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orientations of stretching i.e., material and spatial position of the main axes of 
Green and Cauchy deformation tensors, 

G^VD^V-\ C^UD^U-^, C"^ ^UD-'^U~^. (177) 

Having rigid bodies we can introduce their angular velocities, both "spatial" 
and " CO- moving" versions. Obviously the spatial representations are given by : 

<™' 

respectively for the polar and two-polar decomposition. Similarly, the "co- 
moving" expressions have the known form: 

Obviously, they are related to each other in the usual way, justifying the terms 
spatial and co-moving 

uj,, = Ou,yO-\ Xri = UxriU-\ d,i = vd,iV-^- (180) 

We must remember however that this is something a bit else than the usual 
relationship between, e.g., Q. and Vt. 

The labels "rZ" or "dr" at those angular velocities refer to the fact that those 
quantities describe the L-motion relative with respect to the reference frame 
E, usually non-holonomic one. It is not only geometrically interesting but also 
computationally effective to express the kinetic energy of internal motion Tint, 
cf (|125p through the above quantities. One obtains: 

T.,.. = -lT.(S,7SDVT.(s.7f +iT.(.(f (181) 

where 

W = cDdr = Sdr + O"^— ; (182) 

at 

it is clear that, cDdr is the restriction of 51dr ((52|) to the rigid motion of the 
O-gyroscope, 

Q^B = O-^^ fT'' DcU'' bU^ eV'' . (183) 

Remark: do not confuse completely different things denoted by the same kernel 
symbol S: torsion S^jk, spin S'*j, and deformation-symmetric part of the polar 
decomposition S^^\ unfortunately letters are missing. 

In expression (|18ip we have used the standard orthonormal coordinates in M". 
Then the micromatcrial metric is analytically given by the Kronecker symbol, 

VAB = 5ab- (184) 
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Obviously, if for any reasons convenient, wc can use general rectilinear coordi- 
nates in M" (no (|184[) ). Then in (|18ip instead J we must use J,, given by 

Jri^B^ J'^'^VCB. (185) 

The polar decomposition is convenient when J is general and the internal po- 
tential energy il is spatially isotropic, 

il(L(a;)e(a;)) =il(e(a;)) (186) 

for any isometry L{x) € 0{TxM, g^) at any x G M acting according to (|69ll68p . 
In other words, il depends on e through the Green deformation tensor G (|107p . 
The singular value decomposition (two-polar decomposition p75p Hs convenient 
in doubly isotropic problems, i.e., ones isotropic both in the physical and the 
micromaterial space. This means two things: 

1. Inertial tensor J is invariant under 0{n,rf) acting through (|12[[Til) . there- 
fore, 

jAB ^ j^AB ^* jgAB. (^87) 

the last expression based on the natural choice of orthogonal coordinates 
in M", 

2. Potential energy satisfies both (|186p and 

il(e(a;)i) =il(e(a;)) (188) 

for any L E 0{n,M.) acting through [Ti)) . This means that il depends 
on e only through deformation invariants (jll71 IllSp . 

If those conditions are satisfied, the singular value decomposition provides the 
most effective coordinatization of FM and the formula for kinetic energy be- 
comes: 

Tint = -^Tr (D^x') - ^Tr (t^^^) + /Tr (oxOd) + ^Tr ((^)') , (189) 
where now: 



X = Xdr + Xrl = Xdr + C/ (190) 

Xdr - U-^^fT^dcL^'bL^eV''. 

There is no drive term in i?. Expressions (|1811 1189P have some very peculiar 
features, interesting from the geometric and analytic point of view, and at the 
same time very convenient in physical calculations. Namely, formally they are 
identical with the corresponding formulas for affine motion in flat spaces. The 
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diflFerence is that x ^nd Q contain the additional "drive" terms Xdr, Wdr- These 
terms depend on geometry of M. Translational velocity occurs in these terms, 
because of this, they interfere somehow in (|125p with the translational term Ttr- 
We shall discuss some special examples in two-dimensional spaces n = 2. The 
peculiarity of dimension two is that the group 0(2, R) is Abelian-the exception 
among the groups 0{n, K) because for any n > 2 they are semi-simple. Because 
of this commutativity, the special and co moving representations of angular 
velocity coincide: 

^ Xrl = Xrl, X = X- 

Because of this, another peculiarity appears, namely, there exists an interesting 
class of integrable models with directly separable Hamilton- Jacobi equations. 
However, before doing this, for completeness, we quote the co-moving form of 
balance equations (11481 11491 11521 1153P . They are curved-space afhne counter- 
parts of Euler equations known form rigid body mechanics in Euclidean space. 
Namely, after some calculations one can obtain: 

mh^Bv'' + 2v^v^ScB^ + Is'^dR^'c^bv'' + F^,(191) 

mil''D^''cJ^'^ + N^'', (192) 

where the capitals refer to the co- moving representation, e.g.,w^ = e'*iW*,etc., 
and their raising and lowering is meant in the sense of micromaterial metric tjab 
(usually Sab] rectilinear orthogonal coordinates in M" are most convenient). 
Expressing velocities and affine velocities in terms of linear momentum and spin 
in co-moving representation, we obtain the following balance laws: 

JbcS^^ + -P^'P^'SCB^ + ^^S'^CR^D^BP'' + F^il93) 
m Zm 

E^^JcDS^^'+iV^'', (194) 

The Euler-like structure of those equations is easily readable. 

Obviously, there is nothing wrong in that on the left- hand side of equations we 

have the usual differentiation. Though from the point of view of geometry of 

M, the capital-indices-quantities, i.e., co-moving components, are just scalars. 

Their covariant derivatives along curves are therefore identical with the usual 

derivatives. 



m—— 

dt 



dt 



dp^ 

~dr 

dt 



4 Examples. Two— dimensional homogeneously 
deformable body. 

Let us present as interesting examples the two-dimensional body moving in 
constant-curvature spaces, i.e. , the spherical space 5^ (0, R) and pseudo-spherical 
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Lobachcvsky space ^^'^'+(0, R). There are reaHstic situations when the model is 
rigorously solvable. It is a good illustration of how our method of nonholonomic 
frames works practically. We also obtain an interesting class of models inte- 
grablc in in the Liouville sense. In "polar" coordinates {r,(p) the corresponding 
metric elements are given respectively by 



ds^ = 



+ i?2 sin2 (^) V, 
dr^+i?^ sinh^ d(p^. 



(195) 



In the spherical case all situations with r — 0, r = ttR and arbitrary values of 
(fi correspond to the same points, the "north" pole if r = and the "south" 
pole if r = 2tt. The range of r is [0, nR]. For the pseudo-spherical Lobachcvsky 
space i/^'^'+(0, i?) the range of r is [0,oo]. The most convenient choice of the 
auxiliary reference frame is: 



(r) 



E, 



(r) 



d_ 

dr ' 
d_ 
dr ' 



E, 



iv) 



E, 



1 d_ 
R sin (■^) ' 
1 d 

R sinh (-^r) 



(196) 
(197) 



respectively, in the spherical and pseudospherical case. These frames are evi- 
dently non-holonomic. The kinetic energy may be written in the form: 



(198) 



where are six generalized coordinates. Wc introduce them below; obviously, 
{if, r) will be two of them. Now we can introduce the canonical formalism: 
H = T + V , where T = -^G^^piPj. The matrix G^'' is reciprocal to Gij, i.e. 



where T = ^G^'p.pj 
GijG^'^ = Si''. Translational kinetic energies have the form: 



sphere : Ttr 
pseudosphere : Ttr 




^1 -«'-^'(^) ^ 




.(200) 



When the body is materially isotropic, in the two-dimensional case the inertia 
tensor has only one essential component Ji = J2 = J. It is convenient to use 
the two-polar decomposition (f = UDV~^, where U, V are orthogonal and D-is 
diagonal. Generalized internal coordinates are a, /3, A, /U,: 



U- 



cos (a) — sin (a) 
sin (a) cos (a) 



D 



X 
fj, 



v = 



COS (/?) — sin (/3) 
sin (/3) cos (/?) 
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Then Tint is given by: 

Tint = 



dt) ['dt 



(201) 



where, i? = ^ and x is given by 



da 

X = — + cos 



dt 



on the 2-dimensional sphere and 



da ^ / r \ dtf 



/ r \ d(p 
\r) 'dt 



r\d(fi 



(202) 



(203) 



on the 2-diniensional pseudosphcrc. It is convenient to introduce new variables 
7 := a + /?, 5 := a — P, x := (A — /u), y := (A + /u). Then in the spherical 

case T = ^Gij ^ ^ , where [. . . q'* . . . ] = [r, 7, 6, x, y] and the matrix [Gij] 
is given by 

" 1 " 



[Gij] = 







Gki 







:^ 



and the 3x3 block 



Gki 




is given by 



(204) 



Gk. 



j3:'cos(^) 



Jy^cos(-^) 



In the canonical formalism: T = ^G^^pipj, H = T — V, where GijG^'' = Si''. 
One can show that if generalized coordinates are ordered as [r, (p, 7, 6, x, y], then 
the matrix [G"''] is given by 



[G 



a 61 










(205) 
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where 



fl2 sin2 f 



_cos(j^ 



sir 



TP' J?2sin2(^) 



In the pseudospherical case the analogous equations are vahd. The only differ- 
ence is that trigonometric functions are replaced by (properly signed) hyperbolic 
functions. So we have T = Y^y l^lw"' where [. . .g^ . .] = [r,(p,j,S,x,y] and 
the matrices [Gy] and [G°^] are given in the same form as in spherical prob- 
lem (|204p . (I205p . We must only replace function sin by sinh and cos by 
cosh . In the canonical formalism it is similarly. The generalized coordinates 
are ordered as [r,ip,j,S,x,y], r is now "pseudospherical radius" and r S (O...00). 
Explicitly the kinetic term in the Hamiltonian on the sphere has the form 



T 



Pr^ P% - 2 cos [V^+Ps] 

2m 2TOi?2 sin^ ^ 

(iif^sin2i^+cos2^) {p^+psY 



2J 



Py_ 

2J 



2rni?2 sin 



2 r 



(206) 



2Jj;2 2J?/2' 



and on the pseudosphere 



T 



2 m 



2 cosh (p^ -I- pa) 



2mR? sinh 



cosh^ i ± ^ sinh^ ^ (p, 



W) 



2mi?2 sinh^ ^ 

2 2 
P7 , Ps 



(207) 



2J 2J 2Jx2 2Jy2 



the substitution: S 



-Et 



The Hamilton- Jacobi equation: ^ + H (^q, =0 will be again reduced by 

5*0 (q) ■ The time-independent Hamilton-Jacobi 

E, and we seek for solutions of the form 

So = Sr (r) + (cp) + (7) + Ss (S) + S,{x) + Sy (y) = Sr {r) + lip + C^j + CsS + 
Sxix) + Sy{y). We can separate these equations if ip, 7 and 5 are cyclic variables. 
It is more convenient if we put = i {pa + pp) and ps — ^ {pa — pp) . Then: 
5-0 = Sr (r) + 5^(^) + 5^(7) + Ss{5) + Sx{x) + Sy{y) - Sr{r) + Ip + C^a + 



equation has the form H (^q, 
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CpP + Sx{x) + Sy{y). In the spherical case the Hamilton- Jacobi equation has 
the form: 

2m V rfr / 2mi?2 sin^ i ^ ^ ^ 
In the pseudospherical case: 

2mR? sinh^ ^ ^ ^ 



2m 


\ dr 






2J ( 


y dx J 






2J ( 





8Jx2 

2 



2 

- + 14(x) (209) 



8Jy 



Now we can calculate the action variables. For the all considered cases a, P and 
5 are cyclic variables and the corresponding actions have the same form: 



J = 6 — 't^^dip= I Idip = 2^1 ^ I = , 

J dip Jq 277 

= I ^^^da = Ca da = 27rC„ ^ C„ = — , (210) 

./ da In 2tt 



I fi.,~t ^ \ i 1 1 I 

H = 



Now we can simplify the Hamilton- Jacobi equation using J^, J„, Jg. After this 
we can calculate Jr-,Jx-, Jy The constants of separation and Cy have the 
same form for the sphere and pseudosphere cases: 
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1. On the sphere we obtain: 



J:r. = 



iyJif cos ^ 



. 2m{E-C.^- Cy - Vr{r)) - ^ ^ , " — ^ 



167r^?/^ 



2. On the pseudosphere the formulas read: 



(213) 



J 11 



2m {E-C,-Cy-Vr{r)) 



( - cosh -^) 
47r2i?2 sinh^ 



'2J(C.-K(-))-^^^^rfx, 
/2J(C,-y,(y))-%=^dy. 



(214) 



Substitutmg (|210p into (|210p . (|214p and expressions for Cx and Cj, one obtains 
the exphcit dependence of E on the action variables. The Hamilton- Jacobi 
equation is separable for some realistic (in the elasticity theory sense) potentials 
of the form ii{q) = Vr{r) + Vx{x) + Vy{y) both on the sphere and pseudosphere. 



4.1 Models of potentials 



4.1.1 Models of potentials in the "deformations" plane 



There exists an interesting class of " universally" separable potentials. The cor- 
responding Hamilton- Jacobi equations are separable in all coordinate systems 
used below, i.e., Cartesian coordinates {x,y) and polar coordinates (i;, e) in 
the {x, y)-plane. Among those potentials there are ones effectively integrable, 
realistic, and applicable in elastic problems concerning internal degrees of free- 
dom. Some of them are also separable in elliptic coordinates in the {x, y)-plane, 
however, they are rather non-useful in elastic problems, so we do not make 
use of elliptic coordinates and corresponding action variables. The mentioned 
potentials have the general form: 



Vix,y) = ^ + ^ + C{x^+y^) 
X y 



(215) 
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where A, B, C are constants. An interesting and applicable subclass is given 
by: 

F is constant. The action variables are then given by 

21 



X 



( Jg + JffY 

16 



ttCx 



-7rC„ 



21 F - 



167r2 



Another class of separable potentials very well suited to nonlinear elastic prob- 
lems is given by expressions: 



sin^ (e) ' 

where (<^, e) are polar coordinates in the (x, y)-plane of deformation invariants: 

a; = csin£, j/ = ^cose, (217) 

V^(<;), K(^) are functions of the one indicated variable, respectively e. Ap- 
propriately choosing these "shape functions" we can just obtain the mentioned 
compatibility with the standard requirements of nonlinear elastic dynamics. In 
spherical case (we mean here spherical geometry of two-dimensional "physical" 
space) when generalized coordinates are ordered as [r, a, /3, e, e] , the metric 
[G*-'] , has the form 



1 














Qkl 














m 














m 

J J 
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and 



G 



kl 



^°t'(-H) 

fl2 



where 



^cot2(^) 



fl=sin=(^) 

^cot2(i) 



g°t'(i) 

fl2 



m 



\2 ■ 



J(c cose) 

The separation constants C^, Cy will be combined into a single one, 

C = — Cx — Cy 



(219) 



(220) 
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and 



C = — 



1 (dS,\^ 1 fdSe'" ^ 



The new constant of separation A appears 

^ 27^^; 87rV(sin2£)+*^^^^^' ^^^^^ 



(j2 + 2cos(2e)J,J^ + J2^ 
2J(^ - VM) - ^ , . , -d^^ (223) 

J. = ^ ^2I{C-V,{<;))-^d<;. (224) 
As 1^(£) we propose: 

K(£) = 7COt2(2£), 

where 7 is some constant, then: 

Je = J (^4^/2J(^ + 7)7r - ^8J77r2 + (J„ - J/3)2 - ^8J77r2 + ( J„ + 7^3)2^ . 
As V^(s^) let us take 



where 7 is some constant,then: 



4.1.2 Models for the potential Vr{r) 

The next important problem is to suggest some physically interesting and com- 
putationally effective models for the "radial" potentials Vr{r) in the physical 
space. One can expect that computationally effective will be potentials some- 
how suited to geometry of the "physical" space M. It is well-known that the 
determinant det[5*-'] is a scalar density of weight —2. This is an important geo- 
metric object, used, e.g., for expressing the two-dimensional "volume" (surface 
area) element. Namely, this element is analytically given by 

d/z(r, (f) drdip = J A.et[gij\ drdip. (225) 

-v/det[5(«l V 
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One can expect that simple expresions built of det[(7'^ ] will be good candidates 
for the radial potentials, both mechanically interesting and computationally 
effective. We assume the general form: 

V;(r) =/(r)det[5^^"]. (226) 

/ being something "simple". Obviously, it is constant that is the simplest; we 
denote it by 

/ = Rh, 

7 being also some constant, and then in spherical case where det[(7,y] — j^.^ si,[:j(_r ) 
we have potential Vr{r) in the form: 

7 



sm (-^) 

The "radial" action variable Jr (|210p is then given by: 



Jr = ^/2m4Tr^R^{E + C) + Jl 

1 

2 
1 



-^/2TO47r2i?2^+(J^ + J„)2 (227) 



With —Cx — Cy = C substituted from ()22ip . Obviously, in the compact spherical 
manifold all translational motion are bounded on the "physical" space, therefore, 
the situation where 7 = 0, i.e., Vr{r) — are admissible for the existence of well- 
defined action variables Jr- then we have 



7^\J^ - Jo\ - ^A'h - Jo\ = IJ^ > I Jal I (228) 



/2m4^2^2(^ - - Cy) + J2 - J^. 
In pseudospherical case 

det[(7i,l = = 

^^'^ i?2sinh2(r) 

and again we put 

fir) = 7. 



7 being constant, and then 

V{r) 



i?2 sinh^d^; 
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This also leads to some explicitly integrable expression for Jr (|214p . However 
negative if motion in the r-variable is to be bounded and the corresponding 
action variable Jr to be finite. 

Jr = -^1m^Tx'^R^{C^ + Cy -E) + Jl 
+ iy^2m47r27+(J^ + J„)2 

- iy^2m47r27+ (J,^ - J„)2. 

Solving the formulas (|227p . (|228p with respect to i?, we express the energy as 
a function of all action variables. For that Cx, Cy must be also expressed by 
action variables. One achieves this by solving (|210p with respect to Cx, Cy and 
substituting the resulting expressions to (|228p . (|228p . 



5 Some additional remarks 



We have mentioned above that there exist some models, interesting at least 
from purely geometric point of view of analytical mechanics in itself. This has 
to do with our models of affinely invariant dynamics as developed in [15] , [IS] , 
[5U] . Namely, we can assume in M only some affine connection structure (M, F) 
without any fixed metric tensor g ox\ M . The role of instantaneous metric tensor 
at the configuration e £ FxM G FM will be played by the Cauchy deformation 
tensor; analytically 

C[e],j =77^56^6^,, (229) 

where rj denotes the micromaterial metric. One can always choose coordinates 
in such a way that tjab = Sab (if V is positively definite, what is physically 
assumed). And further on we postulate kinetic energy in the "usual" form: 

= ^VABV^V^ + lvABn^C^''Dj'''', (231) 

where are co-moving components of translational velocity. After some calcu- 
lations based on Poisson brackets we obtain equations of motion in the following 
balance form: 

m {n^b + nb") + 2mw^"S',6° + Y." ^abv'' + F^ (232) 
T."" [VlJ' + -mv^v^ + N''\ (233) 



m 

Dt 

I 

Dt 
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0, (236) 



where indices are raised and lowered with the help of C as a " metric tensor" . 
In the purely covariant and mixed tensor form, these equations acquire the 
following suggestive shape: 

Dn Dv 

-^=m-^ = 2mv\''Scba + ^^dR^cabv'' + F'', (234) 

= N\-mv'^Vb, (235) 

with the connection concerning the shift of indices. These formulas are very 
similar formally to (jl56p . however we must remember about important differ- 
ence: no metric field g is now defined all over M . The Cauchy "metric" C[e] is 
defined only at the point x — 7r(e) where the object is instantaneously present. 
Although no metric field g is fixed in M, nevertheless we can introduce the 
concept of rigid motion by putting the following conditions on motion of our 
object: 

DC^ 
Dt 

and this is simply equivalent to: 

iV" j = = —CjkC ^ i^^i, (237) 

i.e., the C-skew-symmetry of fl. It is interesting that unlike the usual (/-skew- 
symmetry of f2, i.e., the (7-metrical rigid motion, these constrains are in general 
non-holonomic. 

To finish with let us remind our models of affinely-invariant dynamics (in the 
left- and rigid-sense) of extended afhne bodies in flat spaces. There are natural 
in counterparts manifolds. Here we merely quote the corresponding kinetic 
energies, analogous to ones introduced in [48j, [49j, [50| . 
For the spatially affinc and micromaterially isotropic models we have: 

T = jVABV^v^ + ^-VKLv'^^'^^'Mn^ + ^Tr (n^) + |Tr (o)' 

Similarly, for the spatially metrical and micromaterially afhne models (dis- 
cretization of the Arnold description of fluids) we have: 

= '^9^0V'V' + ^5fei5™"^^'m^^'n + J^' ^ + J ■ 

Obviously, in spite of this redundant way of writing, we remember that Tr (f^^) = 
Tr {^^^ 1 Tr [Vl) = Tr (^flj . Models of this kind were analyzed in flat space mo- 
tion [45] , [49] , [50] . In curved manifolds the problem is much more complicated 
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and will be analysed only in some special cases. The more detailed study is 
postponed to the later papers. 
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